clc;
clear;



pde.domain = [0 1 0 1]
pde.meshfile = "mesh_circle_hole.msh"
h = 1/3
bt = ["P1", "P1b", "P2", "P3", "P4", "P5"]
it = ["1d", "2d", "3d"]

mesh = process_mesh(pde)
% fs = process_fespace(mesh, "P2")

% for i = it
% process_showint(i);
% end

% for i = bt
% mesh = process_mesh(pde);

% fs = process_fespace(mesh, i);

% process_showmesh(mesh, fs, [1 1 1]);

% process_showbasis(i, "x");
% end
